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Using methods of high performance computing, we have found indications that knotlike 
structures appear as stable finite energy solitons in a realistic 3+1 dimensional model. 
We have explicitly simulated the unknot and trefoil configurations, and our results sug- 
gest that all torus knots appear as solitons. Our observations open new theoretical 
possibilities in scenarios where stringlike structures appear, including physics of funda- 
mental interactions and early universe cosmology. In nematic liquid crystals and 3 He 
superfluids such knotted solitons might actually be observed. 
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In 1867 Lord Kelvin |TJ proposed that atoms, which at the time were considered as 
elementary particles, are knotted vortex tubes in ether. For about 20 years his theory 
was taken seriously, and motivated an extensive study of knots. The results obtained at 
the time by Tait || remain a classic contribution to mathematical knot theory . More 
recently the idea that elementary particles can be identified with topologically distinct 
knots has been advanced in particular by Jehle [[§]. 

Today it is commonly accepted that fundamental interactions are described by string- 
like structures ||, with different elementary particles corresponding to the vibrational 
excitations of a fundamental string. Even though there are hints of connections between 
modern string theory and knot theory || , knotlike structures have not yet been of much 
significance. 

There are also a number of other scenarios where knotted structures may enter 
Physics 0, 0. These include models in statistical physics, QCD strings that confine 
quarks inside nuclear particles, cosmic strings that are expected to be responsible for 
early universe structure formation, and approaches to quantum gravity where knots are 
supposed to determine gauge invariant observables. Stringlike vortices appear in type-II 
superconductors where they confine magnetic fields within the cores of vortex-like struc- 
tures, and similar phenomena are also present in superfluid 4 He. Recent experiments 
with nematic liquid crystals || and 3 He superfluids have also revealed interesting 
vortex structures that can be described by theoretical methods which are adopted from 
cosmic string models. Finally, physics of knots is rapidly becoming an important part of 
molecular biology, where entanglement of a DNA chain interferes with vital life processes 
of replication, transcription and recombination |jTU |. 

Thus far the physics of knotlike structures has been investigated sparsely. This is 
largely due to a lack of dynamical principles that enable the construction of stable knots. 
One needs a theoretical model where knots emerge as solitons, i. e. as stable finite energy 
solutions to the pertinent nonlinear field equations. 

The literature on solitons is enormous, and there are several extensive reviews H7H, 



Tl|| , |T2| . Until now the activity has mainly concentrated on 1+1 dimensions with the 
notable exceptions of the 2+1 dimensional vortex and nonlinear cx-model soliton, and 
skyrmeons and 'tHooft-Polyakov monopoles in 3+1 dimensions. These are all pointlike 
configurations, and can not be directly associated with knotlike structures. 

When embedded in three dimensions, a pointlike two dimensional soliton becomes 
a line vortex. For a finite energy its length must be finite which is possible if its core 
forms a knot. In 1975 one of us |TJ| proposed that closed, knotted vortices could be 



constructed in a definite dynamical model. The explicit solution suggested in |T3 



is a 



closed torus-like vortex ring, twisted once around its core before joining the ends which 
ensures stability against shrinking. This closed vortex corresponds to the unknot (see 
figure 4) which is the simplest possible knotlike structure. However, despite numerous 
attempts no such stable configuration has been constructed, neither by explicit analysis 
nor by a numerical investigation. 
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Here we shall report on our work to construct knotlike vortices in the model intro- 
duced in [l3||. By employing numerical algorithms in powerful computers we have been 



able to find strong evidence for the existence of the unknot vortex. In addition we have 
found indications supporting the existence of a trefoil vortex (see figure 5). A trefoil 
is the simplest example of torus-knots, which are obtained by winding around a torus 
in both directions |J. Our results suggest that in fact all torus knots should appear as 
vortex solitons in the model proposed in This model describes the 3+1 dimensional 
dynamics of a three component vector ri(x, r) with unit length, n • n = 1. Such a vector 
field is a typical degree of freedom in the nonlinear a-model, a prototype relativistic 
quantum field theory It also appears as an order parameter in the Heisenberg ferro- 
magnet. A unit vector field is also present in models of nematic liquid crystal where it 
characterizes the average direction of the rod, and in 3 He-A superfluid where it deter- 
mines the spin projection direction for a Cooper pair. Indeed, the model proposed in 
T3| is quite universal, and we expect our results to have a large number of applications. 

In order that n(x) describes a localized stationary knot, it must go to a constant 
vector n(x) — » no at large distances. Consequently n(x) defines a mapping from the 
compactified R 3 ~ £ 3 — > S 2 . Such mappings fall into nontrivial homotopy classes 
tt^S 2 ) ~ Z and can be characterized by the Hopf invariant ||. For this we introduce 
the two-form F = (rfnArfn, n) on the target S 2 . Its preimage F* on S 3 is exact, F* = cL4*, 
and the Hopf invariant Qh is given by 



Q H = I FAA (1) 
The Hamiltonian proposed in [ITS 



is 



H = E 2 +E A = J d 3 x g 2 (d fl ii) 2 + e 2 F 2 (2) 

This is the most general three dimensional Hamiltonian that admits a relativistically 
invariant extension in 3+1 dimensions and involves only terms with no more than four 
derivatives. It can be related to the SU(2) Skyrme model when restricted to a sphere 
S 2 G £77(2), but its topological features are different. In particular, the existence 
of nontrivial knotted vortex solutions in (0) is strongly suggested by the lower bound 



H> c- \Q H \i \U 



The first term E 2 determines the standard nonlinear 0(3) a- model which admits 
static solitons in two dimensions. But in three dimensions a scaling x — ► px reveals that 
stable stationary solutions are possible only if £4 is also present. Indeed, under this 
scaling E 2 — > pE 2 but E4 — > p~ x E^ from which we conclude that in three dimensions 
finite energy solutions obey the virial theorem, 

E 2 = E A (3) 
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Several articles have been devoted for analyzing the general properties of the unknot 
vortex in the model @ [Q. However, to our knowledge there have been no real attempts 
to find an actual solution. This is due to the fact, that even in the simplest case of an 
unknot the separation of variables eliminates only one of the three space coordinates, 
leaving numerical methods as the sole alternative for finding a solution. With the recent, 
enormous progress in supercomputing techniques serious attempts are finally becoming 
realistic. 

The Euler-Lagrange equations for (|2|) determine a highly nonlinear elliptic system, 
and a direct numerical approach appears to be complicated. Instead, we consider the 
parabolic equation 

dfa = 5H{<j>) 

dt <50 a 1 ) 

where a denotes a generic dynamical variable in The t-bounded solutions of @ 
connect the critical points of H by flowing away from an unstable critical point at 
t — > — oo towards a stable critical point at t — > +00. From this we expect that by 
defining a suitable initial configuration at t — 0, for large f > we flow towards a stable 
vortex solution of the original stationary equation. 

We parametrize our unit vector by n = (cos (p sin 8, sin ip sin 9, cos 8) through stere- 
ographic coordinates (p = — arctan(^) and 8 = 2 arctan \/U 2 + V 2 . In terms of these 
variables the Hamiltonian @ becomes 

r An 2 Ifip 2 

H = j<?x (i - J + y2)2 {d»U> + dX) + {1 + U i + vr (WW - d v Ud,Vf (5) 

By a global 50(3) rotation we can select our asymptotic vacuum vector rio so that 
outside of the vortex it points to the positive-z direction. This means that outside of 
the vortex we are near the north pole where #(x) pa 0, while the core corresponds to 
the south pole 0(x c ) = it. The internal structure of a vortex can then be investigated 
by cutting it once with a plane at right angle to its core. This cross sectional plane is 
topologically identical to a sphere S 2 : At the core we have 0(x c ) = % corresponding to 
the south pole and outside of the vortex on the plane we have #(x) ~ corresponding 
to the north pole. Furthermore, y(x) increases (or decreases, depending on orientation) 
by 27r when we go around the core once on the cross sectional plane.[] For [/and V this 
means, that we identify them as local coordinates on the Riemann sphere in a patch 
that contains the north pole U(x) = V(x) = 0. 

In figures 1 and 2 we have drawn 8(x) and <£>(x) respectively for our unknot vortex. 
The general structure is clearly visible in these pictures. 

In the present model there are two symmetries that relate a given vortex solution 
with Hopf invariant Qh to an antivortex solution with equal energy but opposite Hopf 

^^More generally, ip is defined modulo 2irn where n is an integer. Here we only consider the simplest 
case with n = ±1. 
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invariant —Qh- The first symmetry emerges when we change orientation cp — ► — (p on the 
cross sectional plane. The second symmetry interchanges the north pole with the south 
pole on the cross sectional plane, i.e. maps the core to 0(x c ) = and the region outside 
of the vortex to 9(x) ~ tt. In terms of the U(x) and V(x) variables this means that 
on each cross sectional plane the Hamiltonian (H) must be invariant under the inversion 
Z(x) = t/(x) + iV(x) — > Z _1 (x). This implies in particular, that the energy density 
//(x) must be concentrated in a relatively narrow tubular neighborhood around the core 
of the vortex. 

In figure 3 we have depicted the energy density on a cross sectional plane for the 
unknot vortex. The concentration of energy density in a tubular manner around the 
core is clearly visible. 

In order to specify the initial condition in (^) we need a knot configuration which is 
topologically equivalent to the desired vortex solution. We specify this configuration by 
first introducing a parametrization x(A) for the center of the knot in R 3 , and then for 
each A use the Serret-Frenet equations to define local coordinates £/(x; A) and V(x.; A) 
on the cross sectional planes. 

The choice of initial parametrization x(A) for the core introduces a scale which may 
be quite different from the one specified by the coupling constant g in In order to 
enhance convergence, we adopt a simple renormalization procedure by promoting g to 
a time dependent function g(t). We define this time dependence by demanding that at 
each value of t the virial theorem (|) must be obeyed. Since a vortex solution obeys (|3|), 
this means that g(t) flows towards a fixed point value g(t) — > g*. 

We have performed our numerical simulations using version 5.3 of G. Sewell's PDE2D 
finite element algorithm fll5]| . During the early phase of our simulations we have used 
the initial knot configuration to determine the boundary conditions on the finite element 
mesh. However, since the exact boundary conditions for the desired vortex are a priori 
unknown, after several time steps we have decreased the size of the finite element mesh 
and used the pertinent simulated configuration to determine the boundary conditions 
in the new mesh. By starting from a sufficiently large initial mesh, we then eventually 
obtain a submesh with boundary conditions that are close to those of the actual solution. 

For our simulation we have used two Silicon Graphics Power Challenge computers 
with R8000 processors equipped with 1GB resp. 2GB of internal memory. 

For the unknot vortex (Qh = ±1), the equations of motion can be simplified using 
axial symmetry. We select the symmetry axis to coincide with the z-axis in R 3 , and 
introduce cylindrical coordinates r, (j), z. With the Ansatz y?(x) = ip(r, z) + and #(x) = 
9(r, z) the ^-coordinate separates, and we obtain a two dimensional equation for U (r, z) 
and V(r, z). That such a separation of variables is possible follows directly from the 
SO(2) x SO(2) symmetry of the unknot configuration. The ensuing equations are defined 
on the half-plane r > 0, z G (— oo, oo), which at each determines our cross sectional 



4 



plane. Besides the inversion invariance, the equations of motion are now also invariant 
under z — > —z. Since the unknot solution is even in z, this halves the CPU time in our 
simulation. 

In figures 1-4 we describe an unknot vortex which is a result of over 50h of CPU time 
with each time step taking about 9 minutes on a 9600 triangle mesh. The mesh has 
been selected so that it is more dense at the boundaries and near the core of the vortex. 
In our simulation we find impressive convergence, allowing us to increase the time step 
by up to 8 orders of magnitude while keeping the relative change in energy intact. The 
Hopf invariant is also very stable, and we have identified the final configuration on the 
basis that the numerically computed Hopf invariant has a slight local maximum with 
Q H = 0.999996 ... 

The trefoil vortex (Qh = ±3) described in figure 5 is a result of almost 200h in 
CPU time on a 21 3 cubic finite element lattice with tri-cubic Hermite basis functions. 
Consequently the number of nodes is about the same as in the case of unknot, but due 
to a lack of any obvious reflection symmetry each time iteration now takes about 20 
minutes of CPU time. The yellow center in figure 5 corresponds to the core x(A) of 
our initial configuration, which has been determined using the energy concept developed 



by J. Simon |lq| . The points in figure 5 have been evaluated using the piecewise cubic 
polynomial approximation obtained from the finite element algorithm. As in figure 4, 
the picture describes a volume where the energy density inside the vortex essentially 
vanishes, i.e. it can be viewed as an extended core of our trefoil vortex. 

We expect that the nonhomogeneity in figure 5 reflects the underlying lattice struc- 
ture that we have used in our simulation. Indeed, a 21 3 lattice is obviously too rough 
to describe our trefoil solution adequately, but unfortunately we do not have access to a 
computer that would allow us to use an essentially larger lattice. Nevertheless, we have 
found definite numerical stability in the sense that the final configuration in figure 5 has 
remained essentially intact under a large number of iterations. We view this stability as 
a strong evidence that we indeed have convergence towards a trefoil vortex solution. 

In conclusion, we have performed numerical simulations with a high-performance 
computer to investigate knotted vortex solutions in the model introduced in |TB[. By 
investigating the unknot and trefoil vortices, we have found strong evidence that torus 
knots indeed appear as solitons. For the unknot, we have found very impressive con- 
vergence and our results for the trefoil are also quite encouraging. However, since the 
computers we can access do not allow us to effectively study dense three dimensional 
lattices, our simulation of the trefoil is still tentative. 

We expect that our results will have numerous important applications. In particular, 
since the order parameter in nematic liquid crystal and 3 He superfluid involves a unit 
three vector, an experimental investigation of vortices in these materials should reveal 
the existence of stable knotlike structures. We also expect that an extension of our work 
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to spontaneously broken Yang-Mills-Higgs theories where stable knotlike vortices can 
not be excluded by scaling arguments, should have interesting physical implications in 
particular to early universe cosmology. 

We wish to thank J.Hietarinta, T.Karki, A.P.Niemi, K.Palo, J.Pitkaranta, J.Rahola, 
D.Riska, R.Scharein, G.Sewell, O.Tirkkonen and J.Tulkki for helpful discussions and 
valuable suggestions. We are particularly indebted to Sami Virtanen for his help with 
visualization and simulating the initial Ansatz for the trefoil, and Matti Grohn for 
helping us with visualization. We are grateful to the Center for Scientific Computing in 
Espoo, Finland for providing us with an access to their line of Silicon Graphics Power 
Challenge computers. 
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Figure Caption 



Figure 1: A combined surface and contour plot of 9(r, z) for an unknot vortex with 
g* pa 0.24, e = 1 in cylindrical coordinates (r, (p, z) and on a cross sectional plane with 
generic cj). The configuration is z ^ —z symmetric. 

Figure 2: A combined surface and contour plot of (p(r, z) for the unknot vortex. 
The line where ip jumps by tt is clearly visible. The configuration is z —>■ —z symmetric. 

Figure 3: A combined surface and contour plot of energy density H (r, z) for the 
unknot vortex, for comparison as figure 1. The configuration is z — > —z symmetric. 

Figure 4: Three dimensional view of the extended core for the unknot vortex. The 
core is defined as the region where energy density essentially vanishes. 

Figure 5: The same as figure 5 but for a trefoil with g* pa 1.3, e = 1. The 
center denotes the core of the initial configuration, determined using a minimum energy 
principle [0. We thank J. Simon for providing us with the initial parametrization. 
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